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Multiplet Structure of Feshbach Resonances in non-zero partial waves 
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We report a unique feature of magnetic field Feshbach resonances in which atoms collide with non- 
zero orbital angular momentum. P-wave (I = f ) Feshbach resonances are split into two components 
depending on the magnitude of the resonant state's projection of orbital angular momentum onto 
the field axis. This splitting is due to the magnetic dipole-dipole interaction between the atoms and 
it offers a means to tune anisotropic interactions of an ultra-cold gas of atoms. A parameterization 
of the resonance in terms of an effective range expansion is given. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The experimental observation of magnetic field Fesh- 
bach resonances (FRs) offers a means by which to widely 
tune the effective interactions in degenerate quantum 
gases. A Feshbach resonance occurs when a quasibound 
state of two atoms becomes degenerate with the free 
atoms and the interatomic potential either gains or loses 
a bound state. As the quasibound state passes through 
threshold the scattering length can be varied in princi- 
ple from positive to negative infinity. FRs were observed 
in bosons fil y, y, 0, y| , in fermions between distinct 
spin states |g, 0, M, Q , and in a single-component Fermi 
gas 10]. Using FRs to study fermions offers a means to 
explore superfluid phase transitions 0, , three-b ody 
recombination H^l m ean-field interactions 0, 0, 0] 
and molecules U^M, 13 E3 • 

Of special interest is the p-wave FR observed in Ref. 
|10|. which exists in a single-component Fermi gas. Due 
to the Pauli exclusion principle the two-body wave func- 
tion must be anti-symmetric under interchange of two 
fermions, implying that only odd partial waves I can exist 
for identical fermions. For I = 1 the Wigner threshold law 
dictates that the p-wave cross section scales as the tem- 
perature squared. This characteristic behavior ordinar- 
ily suppresses interactions at ultra-cold temperatures pl| . 
However, a resonance can dramatically increase the p- 
wave cross section even at low temperatures. 

In this article we discuss characteristics of p-wave Fes- 
hbach resonances. The first is a sensitive dependence of 
observables on temperature and magnetic field. This de- 
pendence arises from a centrifugal barrier through which 
the wave function must tunnel to access the resonant 
state. Only in a narrow range of magnetic field can the 
continuum wave function be significantly influenced by 
the bound state. 

The second characteristic is a doublet in the resonance 
feature arising from the magnetic dipole-dipole interac- 
tion between the atoms' valence electrons. The p-wave 
FRs experience a non-vanishing dipole-dipole interaction 
in lowest order, in contrast to s-wave FRs. This inter- 
action splits the FR into distinct resonances based on 
their partial wave projection onto the field axis, mi = 



or \mi\ = 1. Splitting of the p-wave resonance offers a 
means to tune the anisotropy of the interaction. Dipole- 
dipole interactions in Bose-Einstein condensates and de- 
generate Fermi gas have been considered due to the nov- 
elty of the resulting strong anisotropic interaction p^ . 
P-wave FRs offer an immediately accessible means to ex- 
plore anisotropic interactions in the many-body physics 
of degenerate gases. 

The observed p-wave FR in 40 K occurs when two 
atoms in the \f, m/) = |9/2, — 7/2) hyperfine state col- 
lide. The joint state of the atom pair will be writ- 
ten \hm fl )\f2m f2 )\lmi) = |9/2, -7/2)|9/2, -7/2)|l,m,) 
where mi can take on the values ±1,0. The calcula- 
tions presented here were performed using Johnson's log- 
derivative propagator method p3 ] in the magnetic field 
dressed hyperfine bas is [23 . The potassium singlet and 
triplet potentials |25ll26| are matched to long range dis- 
persion potentials with C% = 3927 a.u. and are fine tuned 
to yield the scattering lengths a s = 104.0 a.u., a t = 174 
a.u. respectively. With these values we are able to re- 
produce the FRs measured in Ref. ||| 0] . 



II. TEMPERATURE AND MAGNETIC FIELD 
DEPENDENCE 

A p-wave resonance is distinct from an s-wave (/ = 0) 
resonance in that the atoms must overcome a centrifu- 
gal barrier to couple to the bound state. The extreme 
dependence on magnetic field and temperature can be 
understood by considering the cross section as a function 
of energy at several values of magnetic field, as shown 
in Fig. H( a )- The lowest curve, with a magnetic field 
B=190 G, shows typical off-resonance p-wave threshold 
behavior. Once the magnetic field is increased to be close 
to the resonance, the cross section changes significantly, 
and a narrow resonance appears at low energy. The res- 
onance first appears for fields just above B=198.8 G. As 
the magnetic field is increased the resonance broadens 
and moves to higher energy. The p-wave resonance's nar- 
rowness is due the fact that atoms must tunnel through 
a centrifugal barrier before they can interact. 

This narrow resonance structure is in stark contrast to 
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the s-wave FR shown in Fig. [2(b), which occurs between 
the spin states |9/2,— 9/2) and |9/2, — 7/2), as reported 
in Ref. 0. The energy dependence of the s-wave FR 
has a much simpler form than the p-wave FR. At high 
energy the cross section is essentially at the unitarity 
limit, which is shown as the solid line. At lower energy 
the cross section plateaus at a constant value of 4ira 2 , 
where a is the s-wave scattering length. The energy at 
which the cross section plateaus depends on the magnetic 
field. The closer to resonance the magnetic field is tuned, 
the lower the energy at which the cross section plateaus. 

The temperature dependence of p-wave FRs results 
from the strong energy dependence of the cross section. 
For a Maxwell-Boltzmann distribution of the atomic en- 
ergies, the thermally averaged cross section is 



(kry 



■ / a{E)Ee- E / kT dE 
Jo 
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where k is Boltzmann's constant and cr(E) is the energy 
dependent cross section pij . 

Figure [21 shows the thermally averaged cross section 
for the mi — 1 p-wave resonance, Fig. [21 (a), and for 
the s-wave FR, Fig. [21(b). The key features of Figure 
[21 (a) are the sudden rise of the cross section and ther- 
mal broadening, which grows dramatically at high field 
as the temperature increases. The rise comes from the 
sudden appearance of the the narrow resonance at posi- 
tive collision energies as the magnetic field is tuned. This 
rise is not temperature dependent because, regardless of 
temperature, the threshold is first degenerate with the 
bound state at a unique magnetic field. By contrast, the 
high field tail of the resonance is sensitive to temperature 
because once the bound state has passed through thresh- 
old the energy dependent cross section peaks at higher 
energies for higher field values. For a fixed magnetic field 
the high field side of the FR, there is a well defined, nar- 
row resonance at a particular energy (Fig. [I] (a)). If the 
temperature is low, very few atom pairs can access this 
resonance. At higher temperatures more atoms experi- 
ence resonant scattering, increasing (a). 

This characteristic asymmetric profile is not present is 
s-wave FRs. The s-wave FR is shown in Fig. [21 (b) near 
its peak. The only effect of temperature in the elastic 
cross section is to wash out the peak of the resonance 
as the temperature is increased. This behavior follows 
from the relatively structureless energy-dependent cross 
section in Fig. (b). 



III. THE DOUBLET FEATURE 

The valence electrons of ultra-cold alkali atoms interact 
via a magnetic dipole-dipole interaction of the form 



2 3(R-si)(R-s 2 )-si-s 2 
-a — (2) 



where a is the fine structure constant, the spin of the 
valence electron on atom i, R is the interatomic separa- 



tion, and R is the normal vector defining the interatomic 
axis. Another way of writing this interaction that isolates 
the spin and partial wave operators is 



Has = 



«V6 
R 3 



£(- 



(si <£>s 2 )-, 



(3) 



Here C^(0,(p) is a reduced spherical harmonic that de- 
pends on the relative orientation of the atoms, and 
(si <g> S2) 2 L q is the second rank tensor formed from the 
rank-1 spin operators j^- C q acts on the partial wave 
component of the quantum state, \lmi), while the Sj's in 
(si <£> S2)_ q act on the electronic spin state of the atoms. 
Equation J3J leads to an interplay of partial wave and 
spin, which contributes an orientation-dependent energy 
to the Hamiltonian. The matrix element of Eq. |J3| in 
our present basis is [24| 
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This term in the Hamiltonian couples different partial 
waves for I' = I ± 2, and it couples different partial wave 
projections mi for I' = I ± 2 and I = I' ^ 0. For elastic 
s-wave scattering (I = V = 0) equation (0} vanishes by 
symmetry. This term only plays a role in s-wave scat- 
tering for s— >d wave transitions. However for p-wave 
scattering (I = I' = 1) this interaction does not van- 
ish, i.e. (lm'i\Cq |lm;) ^ 0. Furthermore, for elastic 
scattering, q — 0, the interaction depends on mj, since 
(11|C#|11) = -~ and (10|Cg|10) = §. The fact that the 
dipole-dipole interaction does not contribute equally to 
all values of mi means that bound states with different 
mi have different energies. This implies that FRs with 
different values of mi couple to distinct bound states and 
thus have different magnetic field dependences. 

The difference between the mi projections can be un- 
derstood intuitively by considering the dipole-dipole in- 
teraction of the two atoms. For the |9/2,— 7/2) spin 
states case the electronic spins are essentially aligned 
with the field. When two dipoles are aligned head to tail 
they are in an attractive configuration, corresponding to 
R-ii = 1 in Eq. (J2J. Vice- versa when the dipoles are side 
by side they are in a repulsive configuration, R • §j = 0. 

Viewing the motion of the atoms in the resonant state 
as classical, circular orbits, the cases mi — and to/ = 1 
are distinguished as in Figure [3 For m; — 0, in Fig. 
[3 (a), the motion of the atoms is in a plane containing 
the magnetic field. Classically this corresponds to mo- 
tion described by the angle 9, where the magnetic field 
lies in the z direction. The interaction for mi — al- 
ternates between attractive and repulsive as the dipoles 
change between head to tail attraction and side by side 
repulsion. On the other hand, for \mi \ = 1, shown in Fig. 
[3(b), the motion of the atoms is in the plane perpendic- 
ular to the magnetic field. Classically this corresponds 
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to motion described by the angle <f>. This interaction is 
only repulsive, because the dipoles are held in the side- 
by-side configuration. Since the dipole-dipole interaction 
for | to/ | = 1 has only a repulsive influence it forms a 
resonant state with higher energy. 

Figure 01 shows the total elastic cross section as a func- 
tion of field at several temperatures. One can clearly 
see the doublet feature in the cross section at low tem- 
perature. The first peak corresponds to |m/| = 1. The 
doublet cannot be resolved at high temperature because 
the width of the resonance is wider than the splitting. 

The energy shift can be estimated using perturbation 
theory. The energy shift due to the dipole-dipole inter- 
action is given in perturbation theory as 

AE mi=0 = - Q y6(10|C 2 |10)(<l> mo ,| (Sl ® 3 82) ° |$m t)(5) 
A£ m(=1 = - a y6(ll|C 2 |ll)($ mo ;| (Sl y 2) ° |0w).(6) 

Here |$ mo i) is the full multi-channel molecular wave func- 
tion without the magnetic dipole-dipole interaction. This 
is the molecular state that couples to the continuum cre- 
ating the FR. We notice that the perturbation is the 
same for each component, except for the angular coef- 
ficients in Equations (0 and (JSJ). When these equations 
are evaluated, we find that the energy difference in the 
molecules, | AE mi= i — AE m[= o\, is 3.7/iK, which is close 
to the closed coupling calculation result of 4.7/iK. As the 
bound states are brought through threshold, we find that 
their energy difference translates into a peak separation 
of 0.5 G, determined from the closed coupling scattering 
calculations. 

Experimentally we have observed the doublet p-wave 
resonance in an ultra-cold gas of 40 K through inelastic 
collisional effects. A gas of atoms in the 1 9/2, —7/2) state 
of 40 K was prepared at T=0.34 in an optical dipolc 
trap characterized by a radial frequency of iv=430 Hz 
and an axial frequency of v z = l Hz \jo\ Il0|. The gas was 
then held at a magnetic field near resonance for 260 [is. 
The resulting Gaussian size of the trapped gas in the 
axial direction was measured as a function of magnetic 
field. The result of this measurement is shown in Fig. [3] 

The observed heating of the gas in Figure is due 
to inelastic processes that occur at the p-wave FR. The 
result clearly shows the predicted doublet structure. The 
splitting between the two peaks is measured to be 0.47±8 
G, in good agreement with the theory. The dominant 
inelastic processes are 3-body losses [HJ, Um > which lie at 
a slightly lower field than the elastic resonance peak ^(J . 

IV. EFFECTIVE RANGE EXPANSION OF THE 
P-WAVE FR 

To compute many-body properties of degenerate gases 
the s-wave scattering length is often used to mimic the 
essential 2-body physics. Near a FR the scattering length 
diverges and can be represented well by a = Ob g (l — 



B _ B ) where db g is the back ground scattering length, A 
is the width, B is the location of the s-wave resonance. 
Scattering length is defined as a =j™ — tan(<5o) jk where 
So is the s-wave phase shift and k — \/2fiE where /i is 
the reduced mass. 

For p-wave collisions the relevant quantity is the scat- 
tering volume, v =i™o — tan(<5i)//c 3 . A simple form like 
the one for a is inadequate for parameterization of the 
p-wave scattering volume because the p-wave resonance 
has a complicated energy dependence. Fig. Efa) shows 
the p-wave scattering volume as a function of field and 
energy. The curves show that as the energy is increased 
the location and width of the resonance change. 

To adequately parameterize the p-wave phase shift 
across the resonance one must use the second order term 
in the the effective range expansion [2^ . Fig. EJb) shows 
A; 3 cot(<5i) plotted as a function of energy for several mag- 
netic fields. This set of curves can then be fit using ef- 
fective range expansion of the form 

fc 3 cot(£i) = -- + ck 2 . (7) 

v 

Where Si is the p-wave phase shift, v is the scattering 
volume, and c the second coefficient in the expansion, 
analogous to the effective range in the s-wave expansion, 
but with units <2q . Both v and c are functions of mag- 
netic field. Fitting v and c to quadratic functions of B, 
which is adequate for the energy range of E < 10~ 6 K 
and magnetic field range of 195 to 205 gauss, we find 

— !— = 8.68155 x 10~ 5 - 8.29778 x 1Q- J B 

Vmi=0 

+1.97732 x 1Q- 9 B 2 
c mi=0 = -1-64805 + 0.015235 - 3.54471 x 10~ 5 .B 2 

— - — = 7.83424 x 10~ 5 - 7.456621 x 10~ 7 B 

«|m,|=l 

+1.76807 x 1Q- 9 B 2 
C| m ,|=i = -2.36792 + 0.022645 - 5.45051 x 10" 5 S 2 

where B is magnetic field in gauss. These fits for i and 
c accurately reproduce fc 3 cot(#i) to within 3% on the in- 
terval specified. This fit does not include experimental 
uncertainties, rather the fit is designed to reproduce the 
the closed coupling calculation with the optimal scatter- 
ing parameters. 

V. IMPLICATIONS FOR EXPERIMENTS 

Because of the p-wave FR, angular dependence of scat- 
tering is under the experimenter's control. For example 
if magnetic field tuned to 199.0 G in the 40 K p-wave 
FR the dominant interactions in the gas will be perpen- 
dicular to the field axis, which follows from the angu- 
lar distribution corresponding to the spherical harmonic 
Yn(#, </>) leading to enhanced collisions rates away from 
the field axis. Whereas if the field is tuned 0.5 G higher, 
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TABLE I: Predicted p-wave resonances in bosonic Rb and 
a Bose-Fermi mixture of Rb-K. The Rb p-wave FRs gain a 
bound state as the field is increased, whereas the opposite 
is true for 40 K p-wave FR. This is why in Rb the mi = 
resonance is lower in field. 



Species Spin States Magnetic Field 



sa Rb 


|2,-2)|2, 


-1) 


B m[ 


=0 


= 247.3 


± 5 G 








B\ m , 


|=i 


= 248.0 


± 5 G 


Hb Rb H/ Rb 


|2,-2) M |1, 


-1>87 


B m[ 


=0 


= 292.8 


± 30 G 










|=i 


= 292.5 


± 30 G 


H 'Rb m K 


|1, l>|9/2,- 


-9/2) 


B m[ 


=0 


= 540.0 


± 30 G 








B\ m[ 


1=1 


= 540.3 


± 30 G 



at cold temperatures the interaction will be dominated 
by Yio(6,4>), characterized by enhanced collisions along 
the field axis. The angular dependence of the collisions 
also has implications for the inelastic 2-body processes. 
These processes are characterized by two atoms gaining 
a predictable amount of energy governed by hyperfine 
splitting and are redistributed in a well defined angular 
manner. 

P-wave FRs offer a means to experimentally study 
anisotropic interactions in systems other than identical 
fermions. For example, we predict that there are p-wave 
FRs in distinct spin states of bosonic 85 Rb and in the 
Bose-Fermi mixture of 87 Rb- 40 K, shown in Table [I] The 
Rb calculations used potentials that are consistent with 
Ref. 29]. The Rb-K calculations are consistent with 



Ref. 30] . On resonance the p-wave cross section be- 
comes comparable to the background s-wave scattering. 
This means that it could have an equally important role 
in determining the collisional behavior and mean-field in- 
teraction of a thermal gas or condensate. 



VI. CONCLUSION 

We have presented characteristics of p-wave FRs. An 
interesting characteristic is the doublet feature for the 
FR with 1 = 1. The splitting is caused by the dipolc- 
dipole interaction having distinct values depending on 
partial wave projection. This might in turn offer a means 
to study anisotropic interactions in quantum gases. An- 
other feature of the p-wave FR is the asymmetric thermal 
broadening, which arises from the resonant state moving 
away from threshold as the magnetic field is tuned. Gen- 
erally, we expect the broadening to occur for fields which 
the v < 0. Just as a p-wave FR splits into two compo- 
nents an Z-wave FR will split into / + 1 components. 
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FIG. 1: (a) P-wave elastic cross section versus energy for 
|9/2, — 7/2>|9/2, — 7/2)|l, 1) collisions for different magnetic 
field values. For each curve the magnetic field in gauss is indi- 
cated. The lowest curve shows an off-resonance cross section, 
(b) For comparison, the s-wave elastic cross section versus 
energy for |9/2, -9/2)|9/2, — 7/2>|0, 0) collisions for different 
magnetic field values. The s-wave FR peaks at B=201.6 G. 
Compared to the p-wave FR these have little structure. The 
solid line is the unitarity limit. 
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FIG. 2: (a) Thermally averaged cross section for 
|9/2, — 7/2>|9/2, — 7/2)|l, 1) collisions as a function of mag- 
netic field. The striking features of this curve are the 
sudden rise and change in width as the temperature 
is increased. (b) Thermally averaged cross section for 
|9/2, — 9/2>|9/2, — 7/2)|0, 0) collisions as a function of mag- 
netic field. The temperature dependence is only evident at 
the peak where it washes out the maximum value. 
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FIG. 3: Schematic representation of classical dipoles interact- 
ing in different circular orbits. Shown in (a) is an orbit with 
mi = 0, which is in a plane containing the magnetic field. 
Here the dipoles sometimes attract and sometimes repel. In 
(b) is shown an orbit with \mi \ = 1, in a plane perpendicular 
to the magnetic field. Here the atoms predominately repel 
one another. 




FIG. 4: The thermally averaged elastic cross section for the p- 
wave FR, including all partial wave projections mi = — 1, 0, 1. 
At low temperatures, the doublet splitting emerges clearly, 
but it is washed out a higher temperatures due to thermal 
broadening. The lower field resonance has \mi\ = 1 and the 
higher field resonance has m ; = 0. 
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FIG. 5: The p-wave FR observed through heating of the gas, 
clearly showing the doublet feature of the p-wave resonance. 
The cloud started at T=0.34 fiK and then was held at a con- 
stant magnetic field. Inelastic processes at the FR, 3-body 
dominated, heat the cloud resulting in an increase in the mea- 
sured size of the trapped cloud. The curve is only a guide to 
the eye. 
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FIG. 6: (a) The p-wave scattering volume for \mi\ = 1 as 
a function of magnetic field at two different energies. Notice 
that both the location at which the scattering volume diverges 
and the width vary with collision energy, (b) fc 3 cot(<5i) for 
the |m;| = 1 p-wave resonance as a function of energy for 
several different values of magnetic field. 



